Dielectronic recombination of xenonlike tungsten ions 
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Dielectronic recombination (DR) of xenonlike W^''"'' forming W^^^ has been studied experimen- 
tally at a heavy-ion storage-ring. A merged-beams method has been employed for obtaining absolute 
rate coefficients for electron-ion recombination in the collision energy range 0-140 eV. The measured 
rate coefficient is dominated by strong DR resonances even at the lowest experimental energies. At 
plasma temperatures where the fractional abundance of W^"^ is expected to peak in a fusion plasma, 
the experimentally derived plasma recombination rate coefficient is over a factor of 4 larger than 
the theoretically-calculated rate coefficient which is currently used in fusion plasma modeling. The 
largest part of this discrepancy stems most probably from the neglect in the theoretical calculations 
of DR associated with fine-structure excitations of the W^''"''([Kr]4d^'' 4/*) ion core. 

PACS numbers: 34.80.Lx,52.20.Fs 



I. INTRODUCTION 

Atomic spectroscopy and collision processes involv- 
ing tungsten ions currently receive much attention, since 
tungsten is used as a wall material in nuclear fusion re- 
actors 0, [1] ■ Consequently, tungsten ions are expected 
to be prominent impurities in fusion plasmas. Radiation 
from excited tungsten ions leads to substantial plasma 
cooling which has to be well controlled in order to main- 
tain the conditions for nuclear fusion. Thus, a compre- 
hensive knowledge of atomic energy levels and collision 
cross sections is required for a thorough understanding of 
the spatial and temporal evolution of the tungsten charge 
states and emission spectra in fusion plasmas To 
date, only a small fraction of the needed atomic data has 
been derived from experimental measurements and most 
comes from theory Q . Because of the complexity of the 
atomic structure of most tungsten ions, the theoretical 
methods require the use of approximations in order to 
become tractable. The associated uncertainties in the 
calculated cross sections are generally difficult to assess. 

The situation is particularly problematic for electron- 
ion recombination which is an important process govern- 
ing the charge balance in plasmas and which also leads to 
the population of excited states in the recombined ion. 
For tungsten ions, no experimental benchmarks in the 
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form of absolute rate coefficients are yet available. Full 
quantum mechanical calculations of rate coefficients for 
dielectronic recombination (DR) have been carried out 
only for a very few charge states 0, i.e., for Ne-like 
W64+ I, g, Ar-like [T^, and Y-like W^^+ [ul- 

Rate coefficients for radiative recombination (RR) have 
been calculated for selected closed-shell and hydrogen- 
like tun gste n ions employing the relativistic Dirac-Fock 
method [H [13. 

Current plasma modeling 0, uses DR rate coeffi- 
cients from the ADAS data base [IJ, [l5| w hich are based 
on the semi-empirical Burgess formula [l6[ . Since the un- 
certainties of these DR rate coefficients were considered 
to be larger than those of rate coefficients for other rele- 
vant atomic processes, the DR rate coefficients were mul- 
tiplied by temperature independent ad-hoc scaling fac- 
tors in order to bring the model predictions into agree- 
ment with measured emission spectra from a tokamak 
plasma [4]. For ions these scaling factors range 

from 0.26 to 2.25, depending on charge state q. 

To the best of our knowledge no experimental measure- 
ments of absolute electron-ion recombination rate coeffi- 
cients of tungsten ions have been carried out so far. Ex- 
periments with highly charged tungsten ions have hith- 
erto mainly focussed on photon emission spectroscopy in 
the extreme ultraviolet and X-ray spectral ranges from 
tokamak and stellarator plasmas 0, [13, [ll] and electron- 
beam ion-traps (EBITs) 1.9 .30]. Results have been ob- 
tained for a number charge states in the range q=27- 
67. These provide vitally needed information for plasma 
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diagnostics. Additionally, cross sections for electron- 
impact ionization have been measured for g=l-10 [3l| . 
and photoionization measurements using synchrotron ra- 
diation have been reported for q—1-3, and 5 (32|. 

Here, we present the absolute experimental rate coeffi- 
cient of Xe-like W^''^ recombining to form Cs-likc W^^+ . 
For these measurements an electron-ion merged beams 
techniqueJ33j| has been employed at a heavy-ion stor- 
age ring [SJl. Details of the experimental procedures are 
given in section |TT1 The experimental results are pre- 
sented and discussed in section IIIIl and conclusions are 
given in section HVl 



II. EXPERIMENT 
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The experiment was performed using the acceler- 
ator and storage-ring facilities of the Max-Planck- 
Institut fiir Kernphysik (MPIK) in Heidelberg, Germany. 
Negatively-charged tungsten ions were produced from 
tungsten carbide in a cesium sputter ion-source and in- 
jected into a tandem accelerator. Positively charged 
W^^"*" ions were obtained by passing the ion beam twice 
(at different energies) through thin carbon foils and by 
appropriate charge-to-mass selection in a dipole magnet. 
The time-averaged electrical current of 186-^20-1- j^j^g g^f_ 
ter this magnet was typically 100 pA at an ion energy 
of 192 MeV. After multiturn injection of five ion pulses 
into the storage-ring TSR, circulating ion currents of the 
order of a few nA were expected considering the stacking 
performance of previous experiments . The velocity of 
the stored ions on their closed orbit was 4.7% of the speed 
of light. Electron cooling [sl] was applied to reduce the 
momentum spread and diameter of the ion beam. Mea- 
surements started 2 s after injection to allow all excited 
states with expected lifetimes of up to several hundred 
milliseconds to decay to the W^°+ ground configuration 
(see discussion below) . 

The TSR storage ring is equipped with two electron- 
beam arrangements, dubbed 'Cooler' and 'Target' [13], 
which can be used for electron cooling and for electron- 
ion collision studies. For the present experiment the Tar- 
get was used for electron cooling and the Cooler served 
as an electron target for the electron-ion recombination 
measurements. This choice was made since the Cooler 
featured a higher electron density than the Target re- 
sulting in larger product-ion count rates. Recombined 
W^^+ ions were separated from the primary W^^"*" beam 
in the first bending dipole magnet behind the Cooler 
and counted with a single particle detector j38|] with ef- 
fectively 100% efficiency. Data were taken at different 
electron-ion collision energies. These were adjusted by 
setting the Cooler cathode voltage appropriately jSS] . 

Usually absolute merged-beams recombination rate co- 
efficients are derived from the measured recombination 
count rates by an appropriate normalization to the elec- 
tron and ion currents [39l . fioj . In the present experiment 
the ion current was so low that it could not be properly 



FIG. 1. (Color online) Number of recombination counts as 
function of storage time t. The Cooler electron-beam was 
switched off at t = 4 s and no Target beam was present at 
any time. The full line is a fit of exponentials to the data 
points. 



measured. Relative recombination rate coefficients were 
determined by normalizing to a proxy of the ion current. 
Here we used the recombination signal at a fixed non-zero 
electron-ion collision energy i?rof- The reference energy 
Ercf — 131.8 eV was chosen such that the recombina- 
tion count rate at this electron-ion collision energy was 
predominantly due to electron-capture during collisions 
with residual gas particles. The reference measurements 
were interleaved with the measurements at energies Ei, 
{ly — 1,2,3,.. .) such that a sequence of electron-ion col- 
lision energies was stepped through during one measure- 
ment cycle (i.e., Ei - S^f -E2- E^of -E3- E.^i -■■■)■ 
The dwell time at each energy was 20 ms. Only the last 
10 ms were used for data taking in order to allow the 
power supplies to settle to their new set values during 
the first 10 ms. Each measurement cycle comprised in- 
jection and cooling followed by a range of 150 energies E^ 
and was repeated for typically one hour. Adjacent energy 
ranges were chosen such that they mutually overlapped 
by 50%. After normalization to the electron current and 
the proxy ion current at E'icf, the data from the vari- 
ous measurement cycles were combined into a relative 
recombination rate coefficient as function of electron-ion 
collision energy. 

The relative recombination rate-coefficient was put on 
an absolute scale by normalization to the recombination 
rate-coefficient ao at zero eV electron-ion collision en- 
ergy. The latter was separately measured by observing 
the storage lifetime of the ion beam j^. To this end 
the decay of the recombination count rate was monitored 
as a function of storage time (Fig. [T]). The lifetime of 
the stored ion beam in the storage ring is limited by 
collisions with residual gas particles. When an electron 
beam is switched on, the lifetime is reduced even further 
by electron-ion recombination. The electron-ion recom- 
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bination rate coefScient can thus be determined from the 
beam decay rate constants A*^""^ and A^°*^^ measured with 
the Cooler switched on and off, respectively, i.e., 

A(°") - Xi°«) 
ao = p77^ — (1) 

where rie — 9.9 x 10^ cm^"^ is the electron density, 
L — 1.5 m is the length of the electron- ion interaction 
region, and C — 55.4 m is the ring circumference. From 
exponential fits to the decay curves shown in Fig. [T] a 
value of ao = (5.3 ± 0.2) x 10~^ cm^ was obtained. 
Repeated measurements reproduced this value within the 
given fit error. Additional uncertainties arise from the 
electron density measurement and the inaccurate knowl- 
edge of the interaction length. A deconvolution of the 
effects of the toroidal electron beam sections on the mea- 
sured rate coefficient [l^ was not carried out since the 
required knowledge of the recombination rate coefficient 
at higher energies is presently not available. The system- 
atic uncertainty of the absolute rate coefficient scale is 
estimated to be 20% at a 67% confidence level [39j . 

The experimental energy spread is mainly determined 
by the velocity distribution of the Cooler electron beam 
which can be characterized by the longitudinal and trans- 
verse teraperatures ksT^i = 0.15 meV and fcsT^ = 
10 meV [39|. These temperature values have been in- 
ferred from a previous experiment which was car- 
ried out under similar conditions as the present measure- 
ments. With these temperatures the experimental energy 
spread is estimated [44| to be 0.04 eV at an energy of 1 eV 
and 0.5 eV at 140 eV. 

For the present measurements no dedicated effort has 
been made to calibrate the experimental energy scale 
more precisely than straight forwardly resulting from the 
merged-beams setup. The resulting systematic uncer- 
tainty is particularly low at very low electron-ion colli- 
sion energies (45j and increases with increasing energy. A 
conservative estimate [s^ yields systematic uncertainties 
of 0.3 and 0.5 eV at electron-ion collision energies of 10 
and 140 eV, respectively. 



III. RESULTS AND DISCUSSION 
A. Merged-beams recombination rate-coefficient 

Figure [2] shows the measured W^'^^ merged-beams re- 
combination rate coefficient as a function of electron-ion 
collision energy. Most dramatically, the rate coefficient 
at energies at least up to 30 eV is characterized by a 
high level about three orders of magnitude above the RR 
rate coefficient estimated from a hydrogenic calculation. 
The measured rate coefficient decreases approximately 
monotonically from eV up to an electron-ion collision 
energy of about 12 eV. From there on, broad resonance 
features become discernible up to the end of the exper- 
imental energy range. The widths of these features are 




Electron-ion collision energy (eV) 

FIG. 2. (Color online) Measured merged-beams rate coeffi- 
cient for electron-ion recombination of W^"^ ions as function 
of relative collision energy. The short-dashed curve is the cal- 
culated RR rate coefficient using a hydrogenic approximation 
[43j . The inset shows the same data in a log-log representation 
and a finer energy binning emphasizing the rate coefficient at 
very low energies. 

much larger than the experimental energy spread. This 
indicates that the observed structures are most probably 
blends of many individually unresolved DR resonances. 

The calculated RR rate coefficient in Fig. [5] has been 
derived from a semi-classical hydrogenic formula for the 
RR cross section [43| with nuclear charge q = 20 and 
cutoff quantum number rimax = 72. This cutoff takes 
into account that recombined ions in loosely bound, high- 
n W^^''"([Kr]4d^° 4/^ nl) Rydberg states are field ionized 
in the TSR dipole magnet before reaching the detector 
(see for details). Although the resulting merged- 
beams RR rate coefficient must be regarded as a crude 
estimate, it is quite clear that the contribution by RR to 
the measured rate coefficient is negligible below ~100 eV. 

The situation is quite similar to what has been ob- 
served in a single-pass merged-beams experiment with 
isoelectronic xenonlike Au^^"*" ions [i^l where the mea- 
sured rate coefficient also exceeded the RR rate coeffi- 
cient by large factors. In the Au^^"*" experiment, two 
well-separated broad groups of DR resonances were ob- 
served in the energy ranges of 15-40 eV and 65-85 eV, 
i.e., in nearly the same energy ranges where there are 
strong DR resonances in the present measurement. How- 
ever the two Au^'''^ resonance groups do not exhibit any 
substructure, whereas the present W^''"'" DR spectrum 
has several distinct local maxima (Fig. [2]). 

The inset of Fig. [2] enlarges the range of the lowest 
electron-ion collision energies where the most obvious fea- 
ture is the rise of the rate coefficient at energies below 
3 meV. This recombination rate enhancement is an arti- 
fact of the electron-ion merged-beams technique. It has 
been studied in detail previously (47j . However, by its 
very limited energy range (0-2 meV in the present scan) 
this feature does not influence the plasma rate coefficient 
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Excitation energy (eV) 

FIG. 3. a) Calculated excitation energies and lifetimes of 
292 out of all 293 excited fine structure components of the 
W^°+([Kr]4d^°4/**) ground configuration. Tlie calculated 
lifetime of the '^Oi2 state is off scale, b) Calculated density 
of states (DOS) comprising all doubly excited 4/* nl states 
above the W^''+(4/'* Fe) ground state with n < 72. 



at the plasma temperatures where W^""'' is expected to 
be formed in fusion plasmas. 

A detailed level assignment of the observed DR res- 
onances is not possible because of the large number of 
excitation channels involved and because of the inherent 
uncertainties of atomic structure calculations for com- 
plex atomic systems. For the interpretation of the low- 
energy DR spectrum in isoelectronic Au^^"*" [i^ , elaborate 
theoretical approaches [48, 49] yielded some quantitative 
agreement with the experimental data but did not iden- 
tify particular states. 

For the present work atomic-structure calculations 
have been carried out using Cowan's atomic structure 
code [scl]. The W^°+ ground state was calculated to be 
the [Kr]4(i^'' 4/^ '^F^ state in agreement with previous 
theoretical results [5l| . Figure [3^ displays the lifetimes 
and excitation energies of all excited states belonging to 
the [Kr]4d^° 4/^ ground configuration. For a given state 
the lifetime was calculated from the Ml and E2 transition 
rates to all accessible energetically lower states. Transi- 
tions of higher order multipoles and configuration mixing 
could lead to a reduction of the calculated lifetimes, but 
were not considered here. 

The majority of levels in the ground configuration have 
lifetimes less than 0.1 s and we expect most of the ion 
beam to relax towards the ground state. However, from 
the present calculations it is found that there are 8 states 
with lifetimes larger than 0.1s. The longest-lived excited 
states are the '^Oi2, ^-fo, and ^Ls states with lifetimes of 
68000, 8.0 and 1.8 s, and excitation energies of 11.10, 
2.43, and 7.276 eV, respectively. Since the waiting time 
after each injection was only 2 s in the present experi- 
ment, this suggests that not all stored ions were in the 
ground state when the recombination measurements took 



place. 

In contrast, the situation in the single-pass Au^^+ ex- 
periment [iHl was quite different. There, the time be- 
tween ion formation and recombination was only about 
5 fis so that most probably a much larger fraction of ions 
were in excited states when they underwent recombina- 
tion in the electron target. Since DR resonance energies 
are different for different initial states, a large fraction of 
various metastable initial states could possibly explain 
the finding that the Au^^"*" DR spectrum is less struc- 
tured than the W^^"*" spectrum. A more definitive state- 
ment about the population of metastable states in the ion 
beams would require an elaborate modeling of the radi- 
ation cascades involved in the de-excitation of the hun- 
dreds of excited states produced by the charge-stripping 
beam-foil interaction. Such a complex theoretical task is 
beyond the scope of the present paper. 

Concluding from the statistical arguments above, the 
presently measured W^°+ DR resonance structure should 
for the most part be associated with excitations of the 
[Kr]4ci^°4/^ ^Fg ground state. Fine structure excita- 
tions to excited states of the ground configuration have 
excitation energies E^xc of up to about 52 eV (Fig. . 
Resonance energies of associated DR resonances can be 
estimated from the Rydberg formula for hydrogenic elec- 
tron binding energies 

E{n) = E,xc - 7^\ (2) 

with TZ 13.606 eV and q = 20, which is reasonably ac- 
curate for sufficiently large principal quantum numbers 
n. Using iSexc = 52 eV one finds that n = 11 is the small- 
est quantum number yielding a positive resonance energy. 
This is well below the cutoff quantum number nmax = 72. 
The smallest excitation energy of the W^"+ ground state 
which yields positive E{n) for n = 72 is 1.05 eV, i.e., 
most of the 293 excited states of the W^""*" ground con- 
figuration contribute to the measured DR spectrum be- 
low 52 eV. The density of states (DOS) comprising all 
(4/^) J nl states with n < 72 which are available for DR 
associated with fine structure excitations within the 4/® 
shell can be easily calculated from the resonance ener- 
gies (Eq. [2]) and statistical weights, i.e. 2rP'{2j + 1), of 
the (4/®)j- nl Rydberg manifolds of DR resonances. As 
can be seen from Fig[3]D the resulting DOS is enormous, 
amounting to more than 10*" states per eV for energies in 
the 5-28 eV range. 

Strong contributions to the observed DR resonance 
structure at energies above 52 eV can be expected from 
the dipole allowed core excitations 4^^^^ 4/* AdPAf^, 
id^Hf 4^1° 4/^ 5d, and 4^1° 4/^ ^ 5g. 
According to the present single-configuration atomic- 
structure calculations the corresponding excitation en- 
ergies are in the ranges 207-238 eV, 157-185 eV, and 
291-318 eV, respectively. Using these excitation energies 
in equation reveals that AcP 4/^ nl resonances with 
n = 5 or 6 and Ad^'^ 4f'^ 5dnl resonances with n = 6 
all have resonance energies in the 0-52 eV range and 
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that resonances with higher Rydberg quantum numbers 
as well as all 46?^" 4/^ 5g nl resonances occur at energies 
above 52 eV. 

It is questionable whether DR resonances associated 
with the just discussed dipole allowed core excitations 
can be held responsible for the steep increase of the ex- 
perimental recombination rate coefficient below 47 eV 
where it rises by almost 3 orders of magnitude from an 
average level of about 2 x 10^^ cm'^ s~^ in the range 
47-80 eV (Fig. [2|). Most probably, a large part of the 
huge low-energy rate coefficient is caused by DR as- 
sociated with fine-structure core excitations within the 
[Kr]4d^° 4/^ ground configuration as represented in Fig. [3] 
and considered in the related discussion. Prominent DR 
resonances of this type have been previously identified 
in DR spectra of much simpler ions such as F-like Fe^^^ 
and 0-like ¥e^^+ 



B. Plasma recombination rate-coefficient 

For the derivation of the W^°+ electron-ion recombi- 
nation rate coefficient in a plasma the measured merged- 
beams rate coefficient (divided by the electron-ion rel- 
ative velocity) has been convoluted with an isotropic 
Maxwellian energy distribution which is characterized by 
the plasma electron temperature T^. In order to obtain 
an accurate result at temperatures far below 1 eV, the 
recombination rate enhancement mentioned above and 
the experimental electron energy spread have to be con- 
sidered carefully [i^. However, these issues do not affect 
the plasma rate coefficient in the temperature range of 
interest for fusion plasmas and, therefore, have been dis- 
regarded here. 

The experimentally derived W^*'^ recombination rate 
in a plasma is shown in Fig. |4l It decreases nearly mono- 
tonically by more than two orders of magnitude across 
the displayed plasma temperature range. At all plasma 
temperatures the experimentally derived rate coefficient 
is very significantly larger than the W^^"*" DR rate coef- 
ficient from the ADAS data base [iBl • The difference 
becomes less drastic with increasing temperature. More- 
over, it should be noted that the present result represents 
only a lower limit because no data were taken at electron- 
ion collision energies higher than 140 eV. Nevertheless, 
the deviation of the experimental result from the theoret- 
ical prediction amounts to a factor of 4.3 at kTc — 160 eV 
where the fractional abundance of W^^"*" is predicted to 
peak in a fusion plasma [s^l (Fig. 2]). 

Interestingly, the agreement between the ADAS 
plasma rate coefficient and the present result becomes 
considerably better at higher temperatures when all DR 
resonances below 47 eV (Fig. [2]) are excluded from the 
convolution procedure that yields the plasma rate co- 
efficient. This finding suggests that the important low- 
energy DR resonances associated with fine-structure core 
excitations have been neglected in the calculation of the 
ADAS rate coefficient. Further support comes from the 
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FIG. 4. (Color online) Left scale: Rate coefficients for the 
recombination of W^"^ ions in a plasma. The thick full curve 
is the present experimentally derived result comprising RR 
and DR with resonance energies below 140 eV. Error bars de- 
note the 20% experimental systematic uncertainty. The long- 
dashed curve results from the experimental data (Fig. [2]) with 
the collision energies restricted to the range 47-140 eV. The 
dash-dotted curve is the recombination rate coefficient from 
the ADAS data base [li,[il. The short-dashed curve is the 
result of a hydrogenic calculation [i^l for the RR rate coeffi- 
cient. Right scale: The thin full line represents the calculated 
fractional abundance of W^"^ in a fusion plasma [sj • 

fact that at low temperatures the ADAS data approach 
the RR rate coefficient (Fig. 01), which has been calculated 
using the hydrogenic cross section mentioned above, but 
with riniax — 1000 instead of rimax = 72. 



IV. SUMMARY AND CONCLUSIONS 

Results from a storage-ring electron-ion recombination 
experiment with a tungsten ion have been presented. Be- 
cause of low ion production rates and high loss rates dur- 
ing electron cooling of stored W^""*" ions, the experimen- 
tal conditions have been less favorable as compared to 
previous experiments with lighter ions. Nevertheless, the 
numbers of mass-selected W^°+ ions stored and cooled 
were sufficient for performing meaningful DR measure- 
ments. The absolute electron-ion recombination rate co- 
efficient, obtained with ±20% systematic uncertainty at 
electron-ion collision energies ranging from to 140 eV. 
An additional error arising from the unknown population 
of metastable states in the ion beam cannot be estimated 
precisely, but does not change our finding that DR via 
fine-structure core excitations is important. 

The experimental recombination spectrum exhibits in- 
dividually unresolved, huge DR resonances at very low 
electron-ion collision energies which strongly influence 
the plasma rate coefficient even at temperatures above 
100 eV, where the fractional abundance of W^""*" is ex- 
pected to peak in a fusion plasma. Because of the ex- 
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traordinary complexity of the W^^"*" atomic structure, no 
definitive assignment of the measured DR resonance fea- 
tures could be made. Atomic-structure calculations sug- 
gest that DR associated with fine-structure excitations of 
the W^'^+([Kr]4(i^° 4/^) ion core makes major contribu- 
tions to the observed low-energy DR resonance strength. 
This fact seems to have been disregarded in the theoret- 
ical calculation of the W^°+ plasma DR rate coefficient 
which is used for plasma modeling by the nuclear fusion 
community. Their resulting rate coefficient is at least a 
factor of 4 lower than the present experimentally derived 
result. 

We also find a factor of 4 difference between our re- 
sults and the ad-hoc scaled W^^"^ DR rate coefficient 
used in recent modeling of a tokamak plasma, for which 
the scaling factor was 0.97 This large discrepancy 

shows that the derivation of rate coefficients from spec- 
troscopic measurements on plasmas generally does not 



yield reliable results. Similar discrepancies should be ex- 
pected for other charge states of tungsten. Consequently, 
the current modeling of the tungsten charge balance in 
fusion plasmas and of the associated radiative cooling 
bears large quantitative uncertainties which can only be 
reduced if further accurate experimental benchmarks for 
the relevant atomic cross sections and rate coefficients in 
the collision energy range below a few hundred eV be- 
come available. 
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